An accurate measurement of electron beam induced displacement cross sections for 

single-layer graphene 
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We present an accurate measurement and a quantitative analysis of electron-beam induced dis- 
placements of carbon atoms in single-layer graphene. We directly measure the atomic displacement 
("knock-on") cross section by counting the lost atoms as a function of the electron beam energy and 
applied dose. Further, we separate knock-on damage (originating from the collision of the beam 
electrons with the nucleus of the target atom) from other radiation damage mechanisms (e.g. ion- 
ization damage or chemical etching) by the comparison of ordinary (^^C) and heavy (^^C) graphene. 
Our analysis shows that a static lattice approximation is not sufficient to describe knock-on damage 
in this material, while a very good agreement between calculated and experimental cross sections is 
obtained if lattice vibrations are taken into account. 
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Radiation damage is one of the key limitations of high- 
resolution transmission electron microscopy (HR-TEM) 
[l|. In particular, the continuous improvements in instru- 
mental resolution [2|, y] inevitably entail increased doses 
per area that need to be applied to a sample. The need 
for high doses is further increased for new techniques 
such as single-atom or single-atomic-column spectroscopy 
y-[7|, atomic resolution electron tomography [S\, or the 
analysis of charge distributions from very high signal-to- 
noise ratio HR-TEM images [9|. For light element mate- 
rials, such as carbon nanotubes lid. Illl. fullerenes Il2l l. 
graphene [ic 



14j . boron nitride 



|16| |. and probably 
many more, the dose limitation is particularly severe for 
three reasons: First, it is obvious, that knock-on dam- 
age cross sections will be higher for low atomic number 
elements [17|. Second, the light elements produce less 
contrast than heavier elements, so that even higher doses 
are needed to obtain a sufhcient signal to noise ratio. 
And third, most of the novel materials from light ele- 
ments, such as graphene or carbon nanotubes, appear in 
the form of low dimensional allotropes that have only one 
or a few atoms in a typical projection of a high-resolution 
image. While almost all atomic spacings can in principle 
be resolved by the currently available instrumentation, 
the question remains whether a sample is stable under 
the beam until an image has been acquired. In spite 
of a wide range of previous studies concerning irradia- 
tion damage in carbon nanostructures |17h24| . a quan- 
titative experimental determination of atomic displace- 
ment cross sections for this important class of materials 
is absent. In fact, only very few quantitative measure- 
ments of electron-beam induced displacement cross sec- 



tions |25l. l2d | (beyond damage threshold measurements 
|27l.l28|) can be found in the literature. The understand- 
ing of irradiation effects is also important for targeted 
irradiation-induced modifications of a material: For the 
case of graphene, for example, a controlled introduction 
of vacancies and non- hexagonal rings may lead to derived 
sp^ hybridized carbon sheets with specific properties [291 - 

Here, we present an extensive measurement and analy- 
sis of electron-beam induced displacements. We directly 
count the number of ejected atoms under irradiation as 
functions of dose, dose rate, and electron energy. Sus- 
pended single layer graphene sheets provide the perfect 
test sample for this analysis: They can be prepared in a 
precisely defined geometry (1 atomic layer thick, hexag- 
onal lattice, with practically no defects initially), are rel- 
atively easy to model, and the number of ejected atoms 
in multi-vacancy configurations can be directly obtained 
from HRTEM images [si,!!^. Under 80 keV electron ir- 
radiation, the defect free graphene lattice remains undis- 
turbed up to very high doses [9, 35] but knock-on damage 
begins already a few keV above this energy |l7H24J . Im- 
portantly, for energies near the knock-on threshold, the 
changes in the lattice occur slowly, so that the appearance 
and growth of multi-vacancies can be directly observed 
in real time. In this way, we can count the number of 
lost atoms as a function of applied dose and for different 
acceleration voltages, hence providing a direct measure- 
ment of the knock-on cross section. 

We present insights from a tremendous data set that 
was obtained for the purpose of quantitating the radi- 
ation damage in graphene: We have obtained and ana- 
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Figure 1: Multi- vacancy defects with increasing dose under 
100 keV observation of ^^C grap Irene (a-c). Images are siiown 
without (left) and with overlay of the atomic configuration 
(right). The dashed line in (a) indicates the area that is used 
to calculate the missing atoms per area. The scale bar is 2 nm. 



lyzed image sequences as shown in Fig. 1 for many accel- 
eration voltages (80, 90, 95, 100 kV), and for both, the 
^^C "normal" graphene sample and isotope-enriched ^"^C 
"heavy graphene" samples. For all of this data, the defect 
configurations were analyzed at different doses of expo- 
sure, and the number of missing atoms was counted (see 
supplementary information for further examples from the 
data set). We also studied ^^C graphene under 20 keV 
electron irradiation [36| , in order to obtain a further dis- 
tinction between knock-on damage and other effects such 
as chemical etching or radiolysis. 

Experimentally, we prepared graphene membranes by 
mechanical exfoliation and transfer to TEM grids as de- 
scribed previously [37|, and by chemical vapor deposi- 
tion (CVD) followed by transfer to TEM grids, as de- 
scribed in Ref. [38]. We assume that these samples con- 
tain the natural isotope composition in carbon, which 
is 98.9% i^c and 1.1% ^^c. In addition, we synthe- 
sized "heavy graphene" samples made from ^'^C, by CVD. 
The synthesis recipe for the ^'^C graphene followed the 
same procedure as we described in Ref. [38|, except that 
the standard methane precursor was replaced by 99% 
^■^C enriched methane (purchased from Sigma- Aldrich). 
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Figure 2: Number of displaced atoms vs. dose and electron 
energy. For the 100 keV case, two different dose rates were 
compared, with (2) having a ca. 3x higher dose rate than 
(1). 



We aligned an image side aberration-corrected FEI Ti- 
tan 80-300 for HRTEM imaging at 80, 90, 95 and 100 
kV, hence providing a closely spaced series around the 
threshold voltage [2Q,]. The spherical aberration was set 
to ca. 20 /im and images were recorded at Scherzer de- 
focus. Under these conditions, dark contrast can be di- 
rectly interpreted in terms of the atomic structure. For 
20 kV imaging, we used an image-side aberration cor- 
rected Zeiss Libra as described in Ref. [36j . In all exper- 
iments, long image sequences of the graphene samples 
were recorded (see supplementary videos), typically con- 
sisting of ~ 100 images with 1 s exposures recorded at 
2-4 s intervals and typical dose rates of 10® '^ ■, . The 
sample is under continuous irradiation, only the beam 
shutter behind the sample is used. We analyze the cre- 
ation and the increase in the density of vacancies and 
multi- vacancies in the image sequences. This approach 
was feasible up to ca. 100 kV, while at 120 kV the damage 
occured too quickly compared to the time or dose needed 
to acquire an HRTEM image with a sufficient signal to 
noise ratio. 

Example images (for 100 kV, ^^C) are shown in Fig. 1.: 
The left hand side shows images from an image sequence 
recorded at 100 kV and the right hand side shows the 
same images with a structure overlay of the atomic con- 
figuration. The analysis of such atomic configurations 
has been described in more detail previously |33l |35|: 
the multi-vacancies reconstruct into configurations that 
involve primarily carbon pentagons, heptagons and oc- 
tagons as well as other non-hexagonal rings, and can 
be well assigned from HRTEM data. For counting the 
atoms, we draw a supercell around the defect clusters, 
such that the boundary of this cell does not intersect any 
defect (see Fig. 1). Moreover, the supercell must not 
contain unpaired dislocation cores, which can be easily 
verified by counting the number of unit cells on opposing 
sides of the parallelogram. We then calculate the num- 



ber of atoms that should be within this cell for the defect 
free case, and compare it to the number of atoms actually 
present. 

The results of this assessment are shown in Fig. 2, 
where the number of lost atoms vs. total dose per area 
is shown for all experiments. A linear fit is made for 
each case, and the slope directly provides the experimen- 
tal knock-on cross section (since a small initial damage 
may be created before the first image is recorded some 
of the lines do not go through the origin). Two inde- 
pendent measurements were made for the lOOkV case, 
but with a 3x different dose rate (3.5 • 10^ -S-^- and 

1 • 10^ nm^s )- From the nearly identical result, we can 
exclude a dose rate effect within our experimental preci- 
sion. Under 80 keV irradiation (not shown in Fig. 2) , no 
vacancies were formed in pristine areas up to very high 
doses (beyond 10^" -^)- 

We begin our discussion by pointing out the clear dif- 
ference between the ^^C and ^'^C graphene membranes, 
and the differences between knock-on damage and a 
chemical etching effect. We find that the generation of 
vacancies within initially pristine, clean and defect free 
graphene membranes depends on the acceleration volt- 
age, and also on the isotope composition (^^C vs. ^■^C). 
Fig. 3a, b shows graphene membranes of the two isotopes 
after exposure to 95 keV electrons, where the difference 
is most clearly visible. Hence, this must be a result of 
a direct collision between a beam electron and the car- 
bon nucleus: Any chemical effect, or ionization damage, 
would not distinguish between ^^C and ^^C. As a side re- 
mark, we note that graphene membranes made from ^^C 
might provide an even better TEM sample support than 
ordinary single-layer graphene: The contrast background 
in HRTEM is identical, but the radiation damage rate is 
lower. 

However, in contrast to the vacancy formation, the 
growth of extended holes in graphene [39j is not predom- 
inantly a knock-on damage effect: We found that the 
growth rate of holes in graphene only weakly depends 
on the electron energy on a wide range of 20 keV to 
100 keV: Holes still form and grow in graphene under 
20 keV irradiation (Fig. Ic-e), and may even grow faster 
at low voltages (see supplementary information). This is 
in stark contrast to expectations from knock-on damage, 
where the threshold for displacing edge atoms is expected 
to be near 50 keV |4fl]. As shown in Fig. [St-d, the ex- 
tended holes always nucleate at contamination sites. We 
noticed that their growth rate is related to the vacuum 
levels, which varied in the range of 10^^ to 10~^ mil- 
libar, e.g. with use of the cold trap in the column, the 
time after insertion of the sample, or different outgassing 
rates of different sample holders. We conclude that this 
is beam-induced etching with residual water or oxygen 
in the system as described in Refs. [23|,|41|. Therefore, 
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Figure 3: Atomic displacements (knock-on damage) vs. chem- 
ical etching. The comparison between ^^C and ^"^C graphene 
shows that the formation of vacancies within the pristine lat- 
tice is a direct knock-on damage effect (a+b, '^^C and ^^C 
sample after a dose of 1.4 • 10^ -^—^ at 95 keV). In contrast, 
the formation of extended holes in graphene is always induced 
by contamination on the sheet (c, initial image, and d after ex- 
posure to ca. 10^ -2-^ ), and the damage rate depends on the 
vacuum levels. Dashed circles in c and d denote same areas. 
(e) Image sequence showing the growth of holes in graphene 
at 20 keV (example shown for a bi-layer area). Scale bars are 
2 nm (a-d) and 5nm (e). 



defect free areas as knock-on effect, but do not take into 
account the extended holes that nucleate at contamina- 
tion sites. 

The analysis of our results culminates in the plot shown 
in Fig. 4. Here, each of the slopes from Fig. 2 provides 
one data point for a measured displacement cross section. 
The error bars indicate the statistical variation (standard 
deviation) in the data. Also shown in Fig. 4 are calcu- 
lated curves from existing and new calculations that will 
be discussed below. For two curves, we show a shaded 
area between Ix and 2x the calculated cross section, 
since correlated sputtering of carbon atoms may increase 
the observed atom loss by up to a factor of two: After 
creation of a mono-vacancy by electron impact, one car- 
bon atom is left with a dangling bond and a much lower 
emission threshold 2l|, |40| . Subsequent (and much more 
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Figure 4: (a) Measured and calculated knock-on displacement 
cross sections. The lower boundary of the shaded areas corre- 
spond to the calculated cross section, while the upper bound- 
ary is twice the calculated value (as would be expected for 
correlated sputtering). Inset shows the calculations for ^^C, 
300K and static lattice on a larger energy range, (b) Corre- 
lated displacement of carbon atoms. After creation of a mono- 
vacancy, one carbon atom remains with a dangling bond and 
a much lower emission threshold. Subsequent sputtering of 
this atom may effectively double the cross section. 



rapid) sputtering of this atom may effectively double the 
rate of atom loss (Fig. 4b) . As a competing mechanism, 
two mono- vacancies that are created close to each other 
may combine and form a stable di-vacancy (since im- 
pacts of energetic electrons can rotate bonds in graphene 
[23, [421 , it is likely that exposure to the beam increases 
also the diffusivity of vacancies). In this case, the sputter- 



ing rate would not have to be doubled. Qualitatively, one 
would expect that correlated sputtering is dominant close 
to the threshold, while nearby mono- vacancies are more 
likely generated at higher electron energies (this is further 
discussed in the supplement). In any case, the result- 
ing multi-vacancy configurations contain only very few 
undercoordinated carbon atoms, while the 3-coordinated 
atoms in the reconstructed configurations are expected 
to have an emission threshold similar to that of an atom 
in the pristine graphene sheet [42]. Another conceivable 
mechanism that might have an influence on the exper- 
imental results, namely the annealing of vacancies with 
mobile carbon adatoms, can not be dominant in our ex- 
periment as evidenced by the absence of a dose rate effect. 
Hence, we expect a rate of atom loss in-between Ix and 
2x of the value calculated for pristine graphene. 

The cross section for Coloumb scattering between an 
electron and a corresponding target nucleus was derived 
by Mott [43|. McKinley and Feshbach have found an 
analytic expression for the Mott scattering cross section 
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where Z is the atomic number of the target atoms. 
En = 13.6ey the Rydberg energy, qq = 5.3 • 10~^^ m 
the Bohr radius of the hydrogen atom, /3 — — (electron 
velocity v^ divided by the speed of light c) , nie the mass 
of th electron and a « r^. T^ax represents the maxi- 
mum transferred energy in the collision event and Tthr a 
threshold energy for atomic displacement. Without mod- 
ifications, equation ^ is suitable to evaluate the total 
"knock-on" cross section for an atom at rest with a given 
ejection threshold energy. The curve from equation ([2]) 
("static lattice" in Fig. 4), features a rather sharp onset 
of radiation damage with increasing acceleration voltage: 
The cross section is zero up to a well defined threshold 
(here 108 kV), and then rises to several barn (beyond all 
our measured values) only a few kV above this thresh- 
old. Changing the displacement threshold in the Mc- 
Kinley Feshbach formula predominantly shifts this curve 
sideways, but does not affect the sharp onset. Hence, 
independent of the free parameter Tthr this approxima- 
tion is in clear contrast to our experiment, which shows 
a smooth onset of the damage cross section between 80 
and 100 keV. 

Remarkably, our data can be explained by considering 
the effect of the struck atom's vibrations on its own dis- 
placement. While the effect has been discussed earlier 
|28l . |45| . |46|, our measurement provides precise experi- 
mental evidence of this intriguing effect. In essence, it 
means that an atom that is struck by an electron while 
it happens to move parallel to the electron beam can 
obtain a higher maximum transferred energy Tmax than 
if it were static. For our calculation, we approximate 
the phonon distribution of the material in the frame- 
work of the Debye model. We use the Debye temper- 
ature calculated for out of plane vibrations in graphene 
of Od = 1287K from Ref. ^. Since 0d depends on the 
speed of sound, it follows for the Debye temperature of 

"C that 0}^ = J^O]^. We extract the distribution of 
atom velocities in the beam direction from the model, 
calculate the maximum transferred energy Tmax{v, E) as 
a function of the atom velocity v and electron energy 
_E, and obtain the weighted sum of the sputtering cross 
section numerically (see supplementary information). In 
other words, we still use the Mott scattering cross sec- 
tion, but we consider that the atom is not at rest initially. 
The threshold energy of Tthr = 22 eV was taken from 
first principles calculations [48] without any adjustments 
(Refs. [2lJ, |22| give similar values). With this value, the 
smooth onset of knock on damage between 80-100 keV is 
very well reproduced. For the first time, no adjustment 
to the calculated threshold energy Tthr is needed to ex- 



plain the data, as was the case in previous studies [22|,|48|: 
Remarkably, the previous mismatch between theory and 
experiments was not due to inadequate calculations of 
Tthr, but because the effects of lattice vibrations on the 
elastic collision were not considered. Interestingly, the 
calculated curves are almost identical for the zero-kelvin 
and room-temperature case (Fig. 4, ^^C T=0 K and 
T=300 K curves). This implies that already the zero- 
point energy of the phonon modes is sufficient to explain 
the increased sputtering cross section as compared to the 
static lattice. 

In summary, we have made an accurate measurement 
of atomic displacement cross sections for carbon atoms in 
single-layer graphene. The cross section smoothly rises 
from practically zero (10^"* barn) at 80 keV to ~ 0.2 barn 
at 100 keV. In practice this means that 80 keV imag- 
ing of defect free graphene is easily possible, while al- 
ready 100 keV TEM images might not represent the orig- 
inal configuration of a sample. A static lattice model is 
not sufficient to model the process, and the contribution 
of atomic motion adds significantly to knock-on damage 
cross sections near the threshold. The difference between 
^^C and ^'^C isotopes is detectable and further confirms 
the model. While the results on graphene will be impor- 
tant for HRTEM studies of this material and related ones 
(especially carbon nanotubes), the generalized insights 
to radiation damage mechanisms should be more gener- 
ally applicable to any material where knock-on damage 
is important. Our results show that knock-on displace- 
ment cross sections can be modeled with high accuracy, 
if lattice vibrations are taken into account. 
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Figures 5-7 show additional image sequences and 
the analysis of defect geometries with the number of 
missing atoms. Figs. 5 and 6 show the radiation damage 
in ^^C graphene under 90 keV and 95 keV electron 
irradiation. For Fig. 6, it should be noted that the 
extended hole that grows on the right-hand side of the 




Figure 5: Radiation damage in ^^C graphene under 90 keV 
electron irradiation. Scale bar Inm. 
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Figurc 6: Radiation damage in ^'^C graphene under 95 keV 
electron irradiation. Scale bar Inm 



image was not taken into account (as discussed above, 
these holes nucleate at contamination sites and grow 
via a beam-induced etching effect). The sequence in 
Fig. 6 also shows two of the rare cases where a heavier 
contamination atom is trapped in a vacancy (these are 
also counted as a lost carbon atom). Fig. 7 shows the 
^■^C graphene sample under 100 keV electron irradiation. 
The rate of atom loss in this sample is significantly 
reduced compared to the ^^C sample (compare Fig. 1 in 
the main article). 



DETAILS OF THE DISPLACEMENT CROSS 
SECTION CALCULATION 

The cross section for Coloumb scattering between an 
electron and a corresponding target nucleus has been de- 
rived by Mott in [43| . McKinley and Feshbach have found 
an analytic expression for the Mott scattering cross sec- 
tion as a function of the maximum transferred energy 
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where Z is the atomic number of the target atoms, i^/j = 
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radius of the hydrogen atom, c the speed of light, (3 = — 
{ve velocity of the electron), c the speed of light, nie the 
mass of the electron, and a = Y37 ■ Tmax represents the 
maximum transferred energy in the collision event and 
Tthr a threshold energy for atomic displacement. For 




Figure 7: Radiation damage in ^''C graphene under 100 keV 
electron irradiation. Scale bar 2 nm 



an atom at rest, the maximuni transferred energy in an 
elastic impact of the beam electron is given by expression 
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where E — ell is the kinetic energy of the electrons 
and rrin the mass of the atomic nucleus. This expres- 
sion is based on the (well justified) approximations that 
the atomic irucleus is much heavier than the electroir 
(to„ ^ TOe), and that its rest energy is much larger than 
the beam energy (c^to„ :s> E). 

The above expression is commonly used to calculate 
sputtering cross sections and thresholds. Our modified 
calculation now explicitly takes into account the fact that 
the target atom is not at rest due to the vibrations of the 
lattice. The importance of lattice vibrations as well as 
thermal effects for the displacement of atoms, irradiated 
with energetic particles was first suggested in 1959 by 
Brown and Augustyniak and subsequently mentioned in 
several publications i28|,|45|,|46|. The effect is not obvious, 
because compared to the energy that is needed to displace 
an atom, its vibrational energy is very small. However, 
the fact that the atom is moving, significantly changes the 
transferred momeirtum and energy upon impact of the 
fast electron. By considering again a (relativistic) elastic 
collision between the beam electron and a moving target 
atom, the maxinuim transferred energy can be written as 
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with: r = \ \/E{E + 2meC^) + rrinV 



and t = ^{E + En){E + 2meC^ + E,,) 
where v is the initial velocity and En ~ Eb^ ig the 
initial kinetic energy of the target atom. Here, v refers 
only to the compoirent of the velocity that is normal to 
the grapheire plane, i.e. parallel to the electron beam. 

We now consider how the cross section (eq. [2]) is 
changed for the case of a moving target. The analyti- 
cal expression for the original Mott Series 1431. 14911 (which 
is accurate up to the medium- Z elements [50|) cair be 
written as [501 
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(5) 
with /3 = -, 6 being the scattering angle and cr^ the 
classical Rutherford scattering cross section: 
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If the target atom is at rest, the angular depeirdence 
of the transferred energy cair be described as 
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(with Tmax as given in eq. [3]) and then integration over 
the scattering angles yields equation [5] 

For a moving target, the angular dependence of the 
transferred energy (with the atom initially moving) can 
be written as: 
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with r and t as given above. This can be rewritten and 
approximated as 

f{e) = f,nax{v,E)*{sm^{-){l-x)+x) « f^a.^sin^ (^) 

(8) 
with X — 1~ -r^TT^j where the approximation is valid for 
a;< 1. 

Hence, we obtain the same expression for ao as in 
equation ^ , except that T^ax (E) is replaced by the new 
Tmax{v, E): For the case of a moving target, we obtain 
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For our parameters (i.e. acceleration voltages, Debye 
temperature, atomic mass) we obtain x w 0.003, meaning 
that the approximation is well justified. The approxima- 
tion can also be understood qualitatively, as that if the 
electron velocity is much higher than the atom velocity, 
the scattering angle is not significantly changed by the 
motion of the atom. Only for very low electron energies 
(on the order of 10 keV or below), one would need to 
consider a modification to equation [2] or [9l 

The mean square velocity of an atom can be calculated 
within the framework of the Debye model to be 
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and it is Gaussian distributed (central limit theorem). 
Here, fcf, is the Boltzmann constant, T the temperature, 
and 6d the Debye temperature. The Debye model natu- 
rally contains a Bose-distribution for the phonon modes 
that also includes the quantum mechanical zero point 
motion for all vibration modes; meaning that the atoms 
are not at rest even at zero temperature. The normalized 
Gaussian distribution yields a probability for finding an 
atom with a velocity between v and v + dv: 

P{v, T)dv = -^=^ exp{^)dv. (11) 



Now, that we have the probability distribution for the 
velocities (fTTjl and using the scattering cross section ((9]) 
with the velocity-dependent expression for T,nax{v) (|4]), 
we can integrate over all velocities weighted by their 
probability in order to obtain the total cross section: 




a{T,E) 



P{v,T)aD{fmaxiv, E))Q{fjnaxiv, E)-Tthr)dv 



(12) 



This last integration step was done numerically. The 
heaviside step function Q{Tjnax{v, E)—Tthr) ensures that 
only positive values for ub are summed up (equations[5]or 
[9]give non-physical, negative values \iTmax < Ttur^ which 
must be replaced by a zero scattering probability). The 

numerical integration was done over a range — 8Vf^ < 
V < Sv w^, providing the curves as shown in Fig. 4 of the 
main article. 



CORRELATED SPUTTERING 

Here, we briefly discuss the two competing mechanisms 
of correlated sputtering and the merging of nearby vacan- 
cies. After removing one carbon atom from the graphene 
lattice, in principle there are three atoms in an under- 
coordinated configuration. However, the configuration 



Figure 8: (a) Cross section for sputtering a 2-coordinated car- 
bon atom compared to 120x the cross section of sputtering a 
3-coordinated carbon (which describes the probability of sput- 
tering a 3-coordinated carbon within a Inm radius around 
the vacancy), (b+c) Radius over which sputtering the un- 
dercoordinated atom would be equally likely as sputtering an 
additional carbon from the pristine lattice, (b) Graphic illus- 
tration for 90 and lOOkV; (c) plot of the radius vs. electron 
energy. 



reconstructs to incorporate a carbon pentagon, which 
closes two of the open bonds and leaves only one atom 
undercoordinated [52]. The displacement threshold for 
this atom is only 14.7 eV [40|, compared to 22 eV for 
an atom in the pristine lattice. Our HRTEM observa- 
tions show predominantly di-vacancies and (clustered) 
multi-divacancies, and exceptionally few single vacancies 
or other undercoordinated atoms. 

Two scenarios appear reasonable in the light of these 
observations. First, one might expect that the under- 
coordinated carbon atom is sputtered shortly after the 
initial vacancy is generated. In this case, the observed 
rate of atom loss would be twice the rate that is calcu- 
lated for the pristine graphene lattice. Second, mono- 
vacancies that form in close proximity to each other may 
'cluster and form di-vacancies without loosing additional 
atoms. As a starting point for quantifying these effects, 
we now compare (a) the probability (cross section) for 
sputtering the undercoordinated carbon atom, and (b) 
the probability for forming two mono-vacancies in close 
proximity. Fig. 8a shows the probabilites (as cross sec- 
tions) for sputtering one undercoordinated atom com- 
pared to sputtering one carbon atom within a Inm ra- 
dius (120 atoms) of the original defect. From the graph 
it is clear that correlated sputtering will be dominant at 
lower electron energies and less significant at high ener- 
gies, but quantification is difficult: Although clustering 
of defects is clearly observed, it is not clear at which rate 
or over which distances it occurs (the choice of Inm as 
the example given above and in Fig. 8a is arbitrary) . Fig 
8b+c shows the analysis represented in a different way: 
It describes a radius over which sputtering an additional 
carbon atom from the pristine lattice would be equally 
likely as sputtering the undercoordinated atom at the 
vacancy. In our data (Fig. 4 of the main article), it in- 



deed appears that the lowest energy data points (-'^^C at 
90kV, ^'^C at lOOkV) correspond to correlated sputtering, 
while ^^C at lOOkV does not. This is in agreement with 
the above discussion if we assume that mono-vacancies 
can cluster into divacancies over a distance of Inm or 
less. In any case, we consider the range from Ix to 2x 
of the calculated cross section in the pristine lattice as a 
(material-specific) uncertainty in the model. 



CHEMICAL ETCHING 

In this section we expand our discussion of chemical 
etching in graphene under TEM imaging conditions. Al- 
though the focus of our analysis is on knock-on damage 
effects, it is also important to identify and separate other 
effects that are not knock-on damage. From a large num- 
ber of experiments on graphene in the TEM, we conclude 
that the growth of extended holes is not predominantly 
an effect of knock on damage. This is particularly sur- 
prising since, under 80 keV irradiation, direct sputtering 
of edge atoms must also occur. However, the experiments 
show that it can not be the dominant mechanism. 

Fig. 9a-c shows the typical formation and evolution 
of a hole in a graphene membrane under 80 keV elec- 
tron irradiation. Initially, a contamination site is present 
(which may be attached to a defect site). Electron- 
energy loss spectroscopy and energy-dispersive x-ray 
spectroscopy indicated that silicon is a predominant con- 
tamination in our samples, and no other elements were 
detected. However, it can not be excluded that other re- 
active species are present in smaller quantities. In any 
case, we almost exclusively observe the formation of holes 
at sites with heavier (than carbon) contamination (Fig. 
9a). The hole rapidly grows while the contamination 
is present and attached to the open edges (Fig. 9b), 
but eventually a hole with pristine (carbon-only) edges 
is formed (Fig. 9c). Such a hole can now be observed over 
extended times and up to high doses [39| . It slowly grows 
under the beam. We have quantified the hole growth 
by measuring the hole perimeter as a function of dose, 
shown for part of our data in Fig. 9d. Hole formation 
and growth was observed in the same way also under 
20kV HRTEM imaging [33]. 

Several observations indicate that the hole growth is a 
chemical reaction with residual gas in the vacuum. The 
most significant one comes from the comparison of 20 
kV and 80 kV data. Knock-on damage of graphene edge 
atoms should be strongly suppressed if not stopped en- 
tirely under 20kV irradiation, compared to 80kV [40|. 
However, we actually measured a slightly more rapid 
growth of holes under the 20kV beam. Second, we found 
that the growth rate of holes strongly varied in differ- 
ent experiments, in spite of similar irradiation conditions. 
Parameters that varied between these experiments were 
the sample holder, the pumping time after mounting the 
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Figure 9: Analysis of chemical etching, (a-c) Typical se- 
quence of hole formation in graphene (scale bar Inm). (a) The 
holes nucleate at a contamination site, (b) The hole rapidly 
grows while heavier contamination atoms are still bound to 
the open edges, (c) Eventually a "clean" hole in graphene is 
formed (only carbon atoms at the edges), and hole growth 
slows down. At this stage, we measure the hole perimeter 
(indicated by the dashed line) and its increase with further ir- 
radiation, (d) Hole perimeter as a function of dose for various 
20kV (red) and 80kV experiments (black). Solid and dashed 
lines denote separate experiments with presumably different 
(although uncontrolled) environmental conditions. 



sample, or the use of the cold trap. In particular, we 
observed that hole growth slowed down after extended 
pumping (several hours after inserting the sample) , com- 
pared to immediately after inserting the sample. We also 
noted one case of a leaking sample holder (contaminated 
o-ring) where hole growth was too rapid for achieving 
any high-resolution images. Indeed, exposure to a low- 
voltage electron beam under a controlled atmosphere of 
water or oxygen was used prev iously for a controlled cut- 
ting of carbon materials |53| . Etching of carbon with 
residual gas in the TEM column was also described by 
other authors [23j. Secondary displacements (e.g. dis- 
placed H atoms that, in turn, displace C atoms) also can 
not play a signifcant role, as again one would expect a 
clear dependence on acceleration voltage for this mech- 
anism. In summary, we conclude that hole formation 
and growth depends on the contamination on the sample 
(which acts as a catalyst) and on the residual vacuum 
levels, and these effects dominate over the knock-on dis- 
placements under our experimental conditions. 



DESCRIPTION OF THE SUPPLEMENTARY 
VIDEOS 



Supplementary video SI: A side- by-side comparison of 
the ^^C and ^'^C graphene sample under 95 keV electron 
irradiation, recorded under identical imaging conditions. 
One can clearly see the much more rapid decay of the 
^^C sample in the electron beam. The total dose at the 
end of the image sequence is 1.7 • 10^-^. 

Supplementary video S2: A long image sequence showing 
a ^^C graphene membrane under 80 keV electron irradi- 
ation. No change can be detected within the pristine 
lattice area. The total dose at the end of the image se- 
quence is 2 • 10^ 
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Supplementary video S3: A long image sequence showing 
a i^C graphene membrane under 80 keV electron irradi- 
ation (showing a larger area view as compared to video 
S2). No change can be detected within the pristine lattice 
area, but the formation of extended holes at the contam- 
ination spots is clearly evident. The total dose at the end 
of the image sequence is 7 • 10^-^. In particular, holes 
start at the heavier (non-carbon) contamination, which 
appears to catalyze the reaction. 

Supplementary video S^: Formation of defects in ^^C 
graphene under 100 keV electron irradiation. The to- 
tal dose at the end of the image sequence is 1.5 • 10^-^. 
A few seconds (less than 30s, total time of the sequence 
is ca. 1 hour) are missing after frames 10, 35, 91 and 132 
as the focus was readjusted. 

Supplementary video S5: Formation of defects in ^"^C 
graphene under 100 keV electron irradiation (note that 
frames 120-177 are out of focus; frames 198-213 are 
recorded under inverse conditions i.e. atoms appear 
white). The total dose at the end of the image sequence 
is 6- 10^-^. 
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